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Abstract 

Modeling,  simulation,  and  integrated 
guidance-control  of  a  kinetic  warhead  utilizing 
moving-mass  actuators  are  discussed.  Moving 
masses  can  be  used  in  any  speed  range  both  in 
the  atmosphere  as  well  as  outside  it,  as  long  as 
there  is  a  force,  either  aerodynamic  or 
propulsive,  acting  on  the  vehicle.  Since  they  are 
contained  entirely  within  the  airframe  geometric 
envelope,  and  because  no  mass  expulsion  is 
involved,  moving-mass  actuation  technique 
offers  significant  advantages  over  conventional 
aerodynamic  control  surfaces  and  reaction 
control  systems. 

The  present  research  developed  a  high 
fidelity,  nine  degree-of-freedom  simulation 
model  of  a  kinetic  warhead  with  three  moving- 
mass  actuators.  This  simulation  model  is  used 
for  actuator  sizing  and  in  the  development  of 
flight  control  systems.  A  software  package  for 
performing  numerical  feedback  linearization 
technique  is  employed  for  the  design  of 
nonlinear  flight  control  systems.  Interception  of 
non-maneuvering  and  weaving  targets  in  both 
atmospheric  and  exo-atmospheric  conditions  are 
demonstrated. 

1.  Introduction 

Some  of  the  earliest  flight  vehicles  were 
controlled  by  moving  the  body  of  the  pilot  [1,2] 
to  effect  a  change  in  the  center  of  mass  of  the 
vehicle.  The  change  in  vehicle  center  of  mass 
alters  the  flight  vehicle  trim,  and  thereby  causes 
a  change  in  the  trajectory.  Advances  in 
aerodynamics  subsequently  made  the  moving- 
mass  approach  to  flight  control  obsolete  in  all 
but  a  few,  specialized  applications.  A  present- 
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day  use  of  moving-mass  control  can  be  found  on 
hang-gliders,  wherein  the  pilot  changes  the 
relative  location  of  his/her  body  with  respect  to 
the  lifting  plane  to  effect  changes  in  flight  path. 
Although  mass  movement  control  is  no  longer 
used  on  commercial  flight  vehicles,  it  continues 
to  be  important  in  applications  such  as  the 
control  of  reentry  vehicles  at  extreme  Mach 
numbers,  where  the  aerodynamic  heating  and 
drag  penalties  can  make  it  impractical  to  deploy 
flight  control  surfaces. 

Recently,  mass  movement  has  been 
proposed  as  a  control  methodology  for  a  kinetic 
warhead  (KW)  in  atmospheric  and  exo- 
atmospheric  engagements  [3].  As  illustrated  in 
Figure  1,  the  moving-mass  control  system 
changes  the  vehicle  center  of  mass  relative  to 
the  external  forces  to  generate  the  desired 
control  moments. 


Figure  1.  Moving-mass  Kinetic  Warhead 
Concept 

For  instance,  if  the  thrust  is  aligned  with  the 
vehicle  longitudinal  body  axis  containing  the 
nominal  center  of  mass  shown  in  Figure  1, 
moving  the  center  of  mass  off  the  body 
centerline  will  result  in  thrust  moments  about 
the  pitch  -  yaw  axes.  Additionally,  roll 
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moments  will  be  generated  if  the  thrust  or  drag 
has  an  angular  misalignment  with  respect  to  the 
longitudinal  axis,  or  if  the  vehicle  is  subject  to 
an  aerodynamic  lift  force.  The  moving-mass 
control  concept  works  equally  well  in  space 
when  the  KW  is  thrusting,  or  in  the  atmosphere, 
when  the  vehicle  experiences  aerodynamic 
forces.  Thus,  this  actuation  technology  can  be 
employed  in  kinetic  warheads  that  have  both 
atmospheric  and  exo-  atmospheric  interception 
capabilities. 

While  the  design  of  flight  control  systems 
using  moving-mass  actuation  appears  to  be 
conceptually  straightforward,  difficulties  arise 
due  to  the  highly  coupled  and  nonlinear  nature 
of  the  system  dynamics.  This  is  partly  due  to 
the  fact  that  in  addition  to  causing  changes  in 
the  vehicle  center  of  mass,  the  moving-mass 
control  system  will  exert  inertial  forces  on  the 
airframe.  Moreover,  the  moving  masses  will 
change  the  instantaneous  moments  of  inertia  of 
the  flight  vehicle,  which  will  then  contribute  to 
changes  in  the  dynamic  response.  The  KW 
control  system  design  must  deliver  the  desired 
interception  accuracy  while  accommodating 
these  dynamic  effects. 

Previous  research  has  demonstrated  that 
integrated  design  of  guidance-autopilot  systems 
can  deliver  fast-responding  flight  control 
systems  by  exploiting  the  synergism  existing 
between  the  vehicle  attitude  dynamics  and  the 
translational  dynamics.  Recently,  integrated 
design  techniques  have  been  investigated  in  the 
context  of  guidance-autopilot-fuze-warhead 
system  design  for  ship-defense  missiles  [4,  5]. 
Three  distinct  methodologies  were  developed 
and  results  were  obtained  for  the  state 
dependent  Riccati  equation  technique  [6  -  8] 
and  the  feedback  linearization  approach  [9  - 
13].  The  research  discussed  in  References  4  and 
5  were  primarily  motivated  by  the  advancements 
in  missile  sensor  and  warhead  technologies.  A 
recent  research  effort  [14,  15]  focused  on  the 
development  of  integrated  guidance-autopilot 
systems  for  a  fixed-aim  warhead  missile. 
Integrated  design  techniques  developed  under 
these  research  efforts  form  the  starting  point  for 
the  present  development  of  moving-mass 
guidance-autopilot  system  for  KW. 


The  objective  of  the  present  research  is  to 
establish  the  feasibility  of  designing  moving- 
mass  actuated  flight  control  systems  for  kinetic 
warheads  that  will  meet  the  accuracy 
requirements  in  both  exo-atmospheric  and 
atmospheric  target  interception  scenarios.  This 
feasibility  demonstration  is  achieved  through  a 
high-fidelity  engagement  simulation  of  a 
moving-mass  controlled  KW  intercepting  a 
tactical  ballistic  missile  in  various  stages  of 
flight. 

Two  different  integrated  guidance-autopilot 
systems  based  on  the  feedback  linearization 
technique  were  synthesized  using  recently 
developed,  computer-aided,  nonlinear  control 
system  design  software  [16].  This  software  was 
developed  under  a  previous  research  effort  [17] 
with  the  Naval  Surface  Warfare  Center, 
Dahlgren,  VA.  A  unique  feature  of  this 
software  package  is  that  it  permits  the  direct 
synthesis  of  nonlinear  control  systems  from 
simulation  models  of  dynamic  systems. 

The  following  sections  will  discuss  vehicle 
modeling,  guidance-autopilot  system  design  and 
evaluation  in  example  engagement  scenarios. 
Detailed  modeling  of  a  KW  with  moving-mass 
actuators  will  be  given  in  Section  2.  Flight 
control  system  designs  and  engagement 
simulation  results  will  be  presented  Section  3. 
Conclusions  from  the  present  research  will  be 
given  in  Section  4. 

2.  Moving  Mass  Actuated  Kinetic  Warhead 
Model 

This  section  will  present  the  formulation  of 
a  9  degree-of-freedom  model  of  the  kinetic 
warhead  with  moving-mass  actuators.  The 
equations  of  motion  were  coded  in  FORTRAN 
and  then  used  in  conjunction  with  Simulink® 
[18]  to  produce  a  simulation  of  the  kinetic 
warhead. 

2.1  Kinetic  Warhead  Model 

In  the  present  study,  the  kinetic  warhead 
consists  of  a  cone-shaped  body  with  three 
actuator  masses  that  can  move  parallel  to  the 
three  orthogonal  axes  of  the  vehicle.  A  free- 
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body  diagram  of  the  kinetic  warhead  is  shown  in 
Figure  2. 

In  addition  to  the  acceleration  due  to 
gravity,  the  kinetic  warhead  will  experience 
aerodynamic  forces  in  atmospheric  flight.  In 
exo-atmospheric  flight,  thrust  is  the  dominant 
external  force  that  acts  on  the  vehicle. 
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are  the  position  vectors  of  the  actuator 
masses  in  the  body  frame  centered  at  the 
vehicle  center  of  mass.  The  mass  offsets 
x 0fjset,  y 0ffset  are  assumed  to  be  specified. 
The  variables  Sx,  5y,  8Z  are  the 
displacements  of  the  three  masses. 
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are  the  relative  velocities  of  the  actuator 
masses  with  respect  to  the  body  B. 
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Figure  2.  Kinetic  Warhead 

2.1.1  Nomenclature 

The  following  notation  is  used  in  deriving 
the  equations  of  motion  (see  Figure  2): 

•  mB  is  the  mass  of  the  kinetic  warhead 
body  B. 

•  mx  mv_  mz  are  the  actuator  masses 
moving  parallel  to  the  body  x,  y,  and  z 
axes,  respectively. 

•  is  the  total  mass  of  the  kinetic 
warhead  and  the  actuator  masses 

•  V0  =  \u  v  w]'  is  the  inertial  velocity 
vector  of  the  center  of  mass  of  the 
vehicle  body  B. 

•  coB  =  [p  q  r]'  is  the  inertial  angular 
velocity  vector  of  the  body  B. 


are  orthonormal  basis  vectors. 

ux,  uy,  u:  are  the  control  forces  acting  on 

the  x,  y,  and  z  axis  actuator  masses, 

respectively. 

IB/o  is  the  3x3  inertia  matrix  of  the 

body  B  about  its  center  of  mass. 

rCP  is  the  position  vector  of  the 

aerodynamic  center  of  pressure  from  the 

center  of  mass  of  body  B. 

rT  is  the  position  vector  of  the  point  of 


application  of  the  external  thrust  vector 
T  from  the  center  of  mass  of  body  B. 

F  is  the  vector  of  aerodynamic  forces  at 
the  center  of  pressure  in  the  body  frame. 
M  is  the  vector  of  aerodynamic 
moments  about  the  center  of  pressure  in 
the  body  frame. 

T  is  the  thrust  vector  in  the  body  frame. 


0 

0 

8 


is  the  gravity  vector  in  the 


inertial  frame.  The  acceleration  due  to 
gravity  is  given  by: 
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f  Re  1 

2 

where  h  is 

mybl2 (  V0-pyxd>B  )+mySy 

kRe  +  h) 

-mybj 

Mb  x  (V0  +0)Bx  py  +  2py )-  CfNBag 

•  =  2.0925x1  ft  is  the  mean 

+  Uy 

g  =  32.174 


radius  of  the  earth. 

•  CNB  is  the  transformation  matrix  from 
the  body  frame  to  the  inertial  frame. 

The  equations  of  motion  of  the  kinetic 
warhead  are  derived  using  Kane’s  method  [19]. 
The  mathematical  model  consists  of  the 
dynamical  and  kinematical  equations  given  in 
the  two  following  subsections. 

2.1.2  Dynamical  Equations 

The  kinetic  warhead  translational  dynamics  are 

given  by  Equation  (1). 

mTV0  -  (mxpx  +  m yPy  +  mzpz )x  cbB 

+  mxpx  +mypy  +mzpz 
=  ~(0B  x  {mTV0  +  mx ( coB  xpx  +  2px  ) 

+  my [o)B  x Py  +2py)+  mz (coB  xpz  +  2pz )} 

+  F  +  nijC^gcig  +  T 

(1) 

The  vehicle  rotational  dynamics  are  given  by 
Equation  (2). 

{mXPx  +  myPy  +  m:Pz)x  VQ  +  lB/0COB 
-  {mxPx  x  Px  +  mypy  X  py  +  mzpz  xp  jx  wB 
+  mxpx  xpx  +  mypy  x  py  +  mzp.  x  pz 
V  -{  „>B  x  I B  0WB  +  mxPx  xaBx  (v0  +aBx  px  +  2px) 

+  rriyPy  x  wB  x  ( VQ  +  coB  x  py  +  2py) 

+  mzp.  xa>Bx  (V0  +  coB  x  pR+  2 pz)  } 

+  M  +(rCPxF)  +  (rmxT) 

+  {mxpx  +  mypy  +  m.p.)x  C‘NBag 

(2) 

The  motions  of  the  actuator  masse; 
described  by  Equations  (3)  -  (5). 


m 


bl(  Vo  -Pxxd,B  )  +  otx4 


x^l 

=  -mxb\ 
+  u  r 


\coB  xfro  +C0B  xPx  +  2Px)-CNBag 

(3) 


(4) 

mzb3  (  V0  ~  Pz  x  cbB  )  +  mjz 
-  ~mzb3  C0B  x  (v0  +  C0B  X  pz  +  2pz)~  C‘NBag 
+  uz 

(5) 

At  each  value  of  the  state  vector,  the  nine 
dynamical  equations  in  (1)  -  (5)  can  be  solved 
as  a  system  of  linear  equations  to  obtain  the 
state  rates  of  the  system  dynamics: 

V0 


coB 

M-A  sx 


-  f(x,u )  where  x  is  the  state 


vector  and  A(x)  is  a  9x9  matrix. 

These  state  rates  are  used  in  conjunction  with 
numerical  integration  algorithms  to  generate 
temporal  histories  of  the  state  variables.  The 
dynamical  equations  have  to  be  augmented  with 
the  kinematical  equations  given  in  the  next 
section  to  complete  the  description  of  the  kinetic 
warhead  dynamics. 

2.1.3  Kinematical  Equations 

The  inertial-frame  position  of  the  mass 
center  of  the  vehicle  body  B  is  given  by 
Equation  (6). 


masses  are 

r  Xs 

l-2( qz  +q23)  2( q,q2  -q3q4)  2(q,q3  +  q2q4  ) 

y 

= 

2(<h<iz  +  q3q4)  !~2(<i2i  +<i23)  itqzcis-qflj) 

(Vo) 

2(<h<l3  ~  ch<i4 )  A  +  q,q4 )  l-  2(q2  +  q2)_ 

~  CNB  Vq 


(6) 


where  q]  -  q4  are  the  Euler  parameters  or 
quaternions. 

The  temporal  evolution  of  the  quaternions  is 
governed  by  the  differential  equation  (7): 
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2.1.4  Kinetic  Warhead  Properties 

The  values  of  the  various  physical 
properties  of  the  kinetic  warhead  are  as  follows. 
These  were  obtained  from  Reference  3 : 

7 lift  =  55  lbm  =  1.7095  slugs 
mx  =  nty  =  mz  =  5  lbm.  =  0. 1554  slugs 
" 0.12821  0  0 

IB/o  =  0  0.6  0  slug-ft^, 

0  0  0.6 


where  the  Ijj  term  was  computed  assuming  a 
solid  right  circular  cone  with  a  base  diameter  d 
=  1  ft,  and  I 22  and  1 33  are  from  Reference  3. 

The  thrust  vector  T  used  for  exo- 
atmospheric  engagements  was  [500  0  0]T  lbf. 


2.1.5  Aerodynamic  Model 


The  normal  and  axial  forces  due  to 
aerodynamics  are  computed  using  the 
expressions: 

FN=QSrefCN,  Fx  ~  ~Q  $  ref  C  a  > 

Q=jpV2,  Sref  -nd2 j4  =  (n/4) 

where  the  coefficients  CN  and  CA  are 
interpolated  from  tabular  data  provided  by  the 
Navy  [20],  and  V  is  the  total  velocity,  the 
magnitude  of  the  inertial  velocity  vector  whose 
body-axis  components  are  u,  v,  and  w.  The 
aerodynamic  coefficients  are  specified  as 
functions  of  Mach  number  and  the  total  angle  of 
attack  aj¬ 


ar  -  cos  —  (9) 

l  V) 

The  axial  and  normal  force  coefficients  plotted 
in  Figure  3  and  Figure  4  illustrate  the  nature  of 
their  dependence  on  Mach  number  and  the  total 
angle  of  attack.  The  normal  force  FN  is  resolved 
into  the  body-axis  y  and  z  components  as: 


Fy  = 


7 


2  2 
V  +w 


fn>  Fz  - 


w 


7 


-F, 


7  7 

v  +w“ 


N 


(10) 

Aerodynamic  moment  coefficients  were  not 
available  at  the  time  of  this  study.  Hence,  the 
aerodynamic  moments  are  computed  using  the 
lift  and  drag  forces  acting  at  the  aerodynamic 
center,  with  the  static  margin  as  the  lever  arm. 
The  static  margin  was  set  to  0.01  ft  for  all  the 
simulations  given  in  this  paper.  The  location  of 
the  aerodynamic  center  is  assumed  to  be 
independent  of  Mach  number  and  total  angle  of 
attack.  Aerodynamic  damping  and  unsteady 
aerodynamic  effects  are  neglected. 


Figure  3.  Axial  Force  Coefficient  vs.  Angle  of 
Attack  and  Mach  Number 


Figure  4.  Normal  Force  Coefficient  vs.  Angle 
of  Attack  and  Mach  Number 


5 


UNCLASSIFIED 


UNCLASSIFIED 


Atmospheric  density  and  the  speed  of  sound 
were  computed  using  the  1976  U.S.  standard 
atmosphere  model.  In  the  present  study,  the 
atmospheric  model  is  defined  for  altitudes  up  to 
435  miles,  although  aerodynamic  effects  are 
expected  to  be  insignificant  above  about  37 
miles. 

2.1.6  Moving-Mass  Positioning  Actuator 
Model 

The  moving  masses  are  assumed  to  be 
positioned  by  proportional  plus  derivative 
servos.  The  servo  gains  are  chosen  to  yield  an 
undamped  natural  frequency  of  50  Hz,  with  0.71 
damping  ratio.  In  response  to  a  commanded 
displacement  of  8ic  the  actuator  servo  will  apply 
a  force 

Fi  -~kp {pi  -8ic)-kd St  (11) 

on  the  moving  mass.  The  proportional  gain  kp 
(lbf/ft)  and  the  derivative  gain  kd  (lbf/ft/s) 
corresponding  to  actuator  natural  frequency  and 
damping,  can  be  computed  as  15340  and  70, 
respectively.  Note  that  these  gains  correspond 
to  a  moving  mass  of  0. 1554  slugs  (5  lb). 

The  displacements  of  the  actuator  masses 
are  limited  since  the  mass  cannot  travel  beyond 
the  confines  of  the  vehicle.  These  constraints 
are  enforced  by  adding  constraint  forces  to  the 
equations  of  motion  so  as  to  bring  the  velocity 
of  the  moving  masses  to  zero  when  the 
displacement  approaches  its  limit  and  to  counter 
the  actuator  force  as  long  as  it  pushes  the  mass 
against  the  stop.  Moreover,  limits  were  also 
imposed  on  the  control  forces  applied  by  the 
servos. 

2.1.7  Line-of-Sight  Angles  and  Rates 

In  addition  to  the  equations  of  motion,  the 
guidance  law  computations  often  require  the 
line -of -sight  angles  and  rates.  The  definition  of 
line-of-sight  (LOS)  angles  Ay,  Az  are  illustrated 
in  Figure  5. 


The  missile  relative-target  position  vector, 
or  the  line-of-sight,  in  the  inertial  frame  is 


r  =  r, 


where 


(12) 

and 

hw  =  \xkw  ykw  zkw J  are  the  inertial 
position  vectors  of  the  target  and  the  kinetic 


t  '  kw 

n  =  h  yt  zi  J 

zkw J  are  the 


warhead,  respectively.  Let  r  -  \rx 


}■ 


The  line-of-sight  angles  can  then  be  defined  as: 


( r  \ 


A}!  —  tan 


r 

\  X  J 


f  „  \ 


,  Az  —  tan " 


(13) 


The  difference  between  the  inertial  velocity 
vectors,  expressed  in  the  Earth-fixed  inertial 
frame,  is  the  relative  velocity,  represented  by 

The  line-of- 


F  m  [o  ry  rzY=Vt-C' 


NBV0  • 

sight  angular  rates  are  obtained  by  taking  the 
time  derivatives  of  the  LOS  angles  above: 


■y  -  rx  tan  Ay 
rxsec  Ay 


■  rz  -  rx  tan  Az 

4  - - j - 


Vy 


~rxry 


rY  sec~  Av 


(14) 


'  X  "z 

J  ,  .. 2 


where  r  =  [rx  +  rF  +  rz  j2 . 

The  expressions  for  angular  acceleration  of 
the  LOS  angle  required  in  one  of  the  integrated 
flight  control  system  designs  can  be  computed 
from  the  following  geometric  relationships. 
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r[rxry  -  rxry )-  2i\rxry  -  rxry  j 

- 7 - ’ 

i  =  (15) 

r3 

where  the  relative  inertial  acceleration  is 
r  =  [rx  ry  r 'z\~at~akw ■  The  inertial 

acceleration  of  the  kinetic  warhead  is 
akw  =  C nb  (To  +  0)  x  Vq),  and  the  acceleration 

of  the  target  a,  is  generally  unknown,  and  needs 
to  be  estimated. 

2.2  Target  Model 

A  non-maneuvering  target  and  a  weaving  target 
model  are  considered  in  the  present  research. 
The  accelerations  of  the  target  in  the  inertial 
frame  are: 

x  -0,  y-0,  z-g  OR  z  —  A  co  ~  sin  a>  t 

(16) 

As  can  be  garnered  from  the  Equation  (16),  the 
nonmaneuvering  target  is  assumed  to  be  falling 
under  the  action  of  gravity.  The  weaving  target 
model  is  assumed  to  be  capable  of  3  g  normal 
acceleration  with  a  frequency  of  0.5  EIz. 

3.  Guidance  and  Control  System 

Design 

The  dynamic  model  of  the  kinetic  warhead 
described  in  Section  2  is  used  to  develop  control 
and  guidance  laws.  This  section  will  present  the 
design  of  guidance  and  control  laws  for 
atmospheric  and  exo-atmospheric  interception 
scenarios.  The  feedback  linearization  technique 
[9  -  13]  is  used  to  design  nonlinear  control  laws 
using  the  complete  model  of  the  kinetic 
warhead.  The  feedback  linearization 
methodology  transforms  a  given  nonlinear 
dynamic  model  into  a  linear,  time-invariant 
form.  Any  linear  control  methodology  can  then 
be  used  to  design  the  control  system.  Inverse 
transformation  of  the  control  law  to  the  original 
coordinates  then  yields  the  nonlinear  control 
law.  The  results  given  in  this  paper  were 
obtained  using  the  pole  placement  design 
technique  in  conjunction  with  the  feedback 


linearization  approach.  A  software  package  for 
feedback  linearization  developed  under  a 
previous  research  effort  with  the  Navy  [17]  is 
used  in  the  design  process. 

The  following  sections  will  discuss  two 
different  integrated  flight  control  system 
designs.  The  integrated  design  methodology 
simultaneously  addresses  the  guidance  and 
autopilot  design  problems.  An  advantage  of  the 
integrated  design  methodology  is  that  it 
eliminates  the  iterations  between  guidance  and 
autopilot  design  processes  required  for 
satisfying  the  flight  control  system  performance 
objectives.  Detailed  discussions  on  the  benefits 
of  the  integrated  design  methodology  are  given 
in  References  4,  5,  14,  15,  21  and  22. 

3.1  Integrated  Design  Based  on  Line-of- 
Sight  Rate  Regulation 

This  integrated  design  approach  is  based  on 
the  proportional  navigation  concept,  in  the  sense 
that  the  flight  control  system  is  designed  with 
the  objective  of  driving  the  line-of-sight  rates  to 
zero.  The  integrated  guidance-autopilot  system 
is  also  responsible  for  stabilizing  the  kinetic 
warhead  dynamics.  In  order  to  design  the 
integrated  flight  control  system  using  the 
Nonlinear  Synthesis  Tools  [16]  software 
package,  the  control  influence  chain  can  be 
defined  as: 

uz  — »  8Z  — »  Sz  —>  q  — »  w  — »  Xz  (17) 

Uy  Sy  8y  >'  V  Ay  (18) 

These  relationships  describe  the 
interdependence  of  control  and  state  variables  in 
the  kinetic  warhead  dynamics.  For  instance,  the 
first  expression  suggests  that  the  force  uz 
acting  on  the  moving  mass  influences  the 
velocity  A,  and  position  A,  of  the  moving  mass 
along  the  z-axis,  which  in  turn  influences  the 
kinetic  warhead  pitch  rate  q  resulting  in  a 
change  in  the  body  velocity  component 
w  along  the  z-axis.  This  velocity  component 
will  result  in  an  acceleration  component  normal 
to  the  kinetic  warhead,  which  will  then 
influence  the  line-of-sight  rate  Xz .  A  similar 
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description  can  be  also  be  generated  for 
Expression  (18). 

The  control  influence  chains  are  provided  to 
the  nonlinear  control  system  design  software  in 
the  form  of  a  matrix,  Reference  16  gives 
additional  details  on  how  this  matrix  can  be 
formulated  in  a  given  problem. 

3.1.1  Engagement  Scenario  1 

The  first  engagement  scenario  used  to 
evaluate  the  integrated  guidance-autopilot 
system  is  described  below.  The  closed-loop 
poles  in  both  channels  were  {-51,  -50,  -35,  -30, 
-20}.  The  kinetic  warhead  and  target  trajectories 
are  shown  in  Figure  6  and  Figure  7. 

This  engagement  scenario  is  similar  to  the 
endo-atmospheric  engagement  considered  in 
Reference  3  in  which  both  the  warhead  and 
target  are  initially  at  an  altitude  of  45,000  ft,  and 
50,000  ft.  apart.  Both  vehicles  have  15-degree 
flight  path  angles  and  are  on  reciprocal 
headings,  and  both  have  an  initial  velocity  of 
6000  ft/sec.  In  addition,  the  target  has  an  initial 
offset  of  1000  ft.  in  the  east  direction . 

The  miss  distance  for  this  engagement  was 
0.28  ft.  The  line  of  sight  rates  are  shown  in 
Figure  8.  The  integrated  guidance-autopilot 
drives  the  line-of-sight  rates  to  near-zero  early 
in  the  engagement,  and  then  attempts  to  keep  it 
there.  Note  that  it  is  possible  to  introduce 
dynamic  compensation  networks  in  the 
integrated  guidance-autopilot  to  distribute  the 
control  effort  more  uniformly  over  the 
engagement.  Moreover,  the  gain  in  the  line-of- 
sight  rate  states  can  be  increased  to  further 
reduce  the  miss  distance. 


Figure  6.  Warhead  and  Target  Trajectories, 
Vertical  Plane 

The  angle  of  attack  and  angle  of  sideslip  are 
shown  in  Figure  9.  There  is  a  large  initial 
acceleration  seen  in  Figure  10,  particularly  in 
the  lateral  channel,  to  remove  the  LOS  error.  A 
noticeable  roll  rate  is  induced  due  to  the  moving 
masses,  which  can  be  seen  in  the  angular 
velocities  shown  in  Figure  11.  The  actuator 
mass  displacements  and  the  control  forces  are 
shown  in  Figure  12. 


Figure  7.  Warhead  and  Target  Trajectories, 
Horizontal  Plane 
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Figure  8.  Line  Of  Sight  Rates  vs.  Time 


Figure  9.  Angles  of  Attack  and  Sideslip  vs. 
Time 


Figure  11.  Body  Angular  Velocities  vs.  Time 


Figure  10.  Body-Frame  Accelerations  vs. 
Time 
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Figure  12.  Actuator  Mass  Displacements  vs. 

Time 

3.1.2  Engagement  Scenario  2 

The  second  scenario  is  the  exo-atmospheric 
case  from  Reference  3  where  the  warhead  and 
target  start  at  an  altitude  of  545,000  ft.  and  are 
120,000  ft.  apart  on  reciprocal  headings  and 
flight  path  angles  of  45  degrees.  Both  vehicles 
have  an  initial  velocity  of  6000  ft/sec.  In  this 
case  the  poles  had  to  be  slowed  down  due  to  the 
smaller  force  available  for  steering.  The  closed 
loop  poles  were  {-51,  -50,  -25,  -20,  -10}  in  both 
channels.  The  kinetic  warhead  and  target 
trajectories  are  shown  in  Figure  13.  The  miss 
distance  was  0.42  ft. 

The  line  of  sight  rates  are  shown  in  Figure 
14.  As  in  the  previous  engagement  scenario,  the 
integrated  guidance-autopilot  system  proceeds 
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to  correct  most  of  the  line-of-sight  rates  at  the 
beginning  of  the  engagement.  The  angle  of 
attack  and  angle  of  sideslip  are  shown  in  Figure 
15.  The  accelerations  are  seen  in  Figure  16. 
The  accelerations  shown  are  in  the  wind  frame, 
i.e.  perpendicular  to  the  velocity  vector,  instead 
of  the  body  frame  as  was  shown  in  the 
atmospheric  case.  This  is  because  of  the  fact 
that  since  the  maneuvering  forces  are  derived 
primarily  from  the  thrust  vector  acting  along  the 
kinetic  warhead  longitudinal  axis,  acceleration 
components  normal  to  the  velocity  vector  can  be 
nonzero,  while  the  lateral  and  normal 
acceleration  components  along  the  body  axes 
remain  near  zero. 

As  in  the  atmospheric  interception  scenario, 
there  is  a  large  increase  in  acceleration  at  the 
beginning  of  the  engagement.  Note  that  since 
the  thrust  is  500  lbf  and  the  total  mass  is  70  lbm, 
the  maximum  acceleration  obtainable  is  around 
7  g’s.  Figure  15  shows  that  the  maximum  angle 
of  attack  during  the  engagement  was  about  68 
degrees.  Note  that  this  fact  has  important 
bearing  on  the  types  of  seekers  that  can  be 
employed  in  conjunction  with  the  moving  mass 
control  of  kinetic  warheads.  Engagements  can 
be  accomplished  with  smaller  angle  of  attack 
and  angle  of  sideslip  if  higher  thrust  levels  are 
available.  The  angular  velocities  and 
quaternions  are  shown  in  Figure  17.  The 
actuator  mass  displacements  and  the  control 
forces  are  shown  in  Figure  18. 


Figure  13.  Warhead  and  Target  Trajectories, 
Vertical  Plane 


Figure  14.  Line  of  Sight  Rates  vs.  Time 


Figure  15.  Angles  of  Attack  and  Sideslip  vs. 
Time 


Figure  16.  Wind-Frame  Accelerations  vs. 
Time 


10 


UNCLASSIFIED 


UNCLASSIFIED 


Time  (sec) 


Figure  17.  Body  Angular  Velocities  vs.  Time 


Figure  18.  Actuator  Mass  Displacements  vs. 
Time 


3.1.3  Engagement  Scenario  3 

This  engagement  features  a  weaving  target 
example  from  Reference  3.  The  warhead  and 
target  are  initially  at  an  altitude  of  45,000  ft.  and 
are  headed  directly  toward  each  other.  The 
initial  velocities  are  6000  ft/sec.  for  both 
vehicles.  The  target  weaves  in  the  vertical  plane 
with  an  acceleration  of  3  g  and  a  frequency  of 
0.5  Hz.  The  poles  of  the  integrated  flight 
control  system  were  chosen  to  be  the  same  as  in 
Subsection  3.2.1. 

The  kinetic  warhead  and  target  trajectories 
are  shown  in  Figure  19.  The  terminal  miss 
distance  for  this  case  was  0.28  ft.  The  line  of 
sight  rates  are  shown  in  Figure  20.  As  in  the 
previous  engagement  scenarios,  the  flight 


control  system  attempts  to  maintain  the  line-of- 
sight  rates  close  to  zero.  However,  the  flight 
control  system  encounters  difficulties  in 
maintaining  these  quantities  close  to  zero  due  to 
the  dynamic  nature  of  the  target.  The  angle  of 
attack  and  angle  of  sideslip  are  shown  in  Figure 
21.  The  actuator  mass  displacements  and  the 
control  forces  are  shown  in  Figure  22. 


Figure  19.  Warhead  and  Target  Trajectories, 
Vertical  Plane 


Figure  20.  Line  of  Sight  Rates  vs.  Time 
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Figure  21.  Angles  of  Attack  and  Sideslip  vs. 
Time 


Figure  22.  Actuator  Mass  Displacements  vs. 

Time 

3.2  Integrated  Design  Based  on  Collision 
Course  Tracking 

The  second  version  of  the  integrated  flight 
control  system  design  continuously  attempts  to 
align  the  kinetic  warhead  velocity  vector  along 
the  collision  course  with  the  target.  This  is 
accomplished  by  constructing  a  tracking  system 
for  the  flight  path  and  heading  angle  commands 
and  a  command  generator  that  computes  the 
necessary  angles  to  achieve  a  collision.  The 
flight  path  angle  (y)  and  the  heading  angle  (%) 
are  defined  as: 


y  -  -  tan 


-1 


tJx 


2  .  ■ 2 


x'  +  y 


,  x  = tan 


-1 


fZ 

U' 


(19) 

where  [x  y  z  ]*■  =  CjygVQ  are  the  kinetic 

warhead  velocity  components  expressed  in  the 
inertial  frame.  The  errors  between  the  actual 
and  desired  angles  are  integrated  with  time,  and 
these  integrator  states  are  also  fed  back.  The 
control  influence  chains  for  the  integrated 
control  system  design  are: 


uz^Sz^Sz^q^w^y^ \\yc  -  y]dt 

(20) 

Uy  -►  Sy  ->  Sy  T  V  J  fac  ~  Z^ 


(21) 

The  collision  course  is  computed  by 
assuming  that  the  target  and  warhead  fly  at 
constant  speeds  and  directions  until  collision  at 
time  tg0.  The  kinetic  warhead  velocity 

components  required  to  achieve  interception  can 
be  computed  from  the  collision  conditions.  For 
instance,  along  the  x -direction  in  the  inertial 
frame,  the  collision  conditions  require: 


HW  X'  /  ■  gl)  KW  ”KW 

—  Xf  {t go ) —  (p)  +  Igo^t 


(22) 


This  expression  can  be  simplified  to  yield: 
•  _xt(0)-xkw(0) 

xkw  ,  '  xt 


lgo 


(23) 


The  kinetic  warhead  velocity  vector  components 
along  the  y  and  z  directions  can  be  computed 
similarly.  If  information  about  the  target 
acceleration  is  available  it  can  be  included  in  the 
above  equation.  The  time-to-go  is  estimated 
from  the  range  and  range  rate: 

lgo  -  (24) 

r 

The  differential  equations  for  flight  path  and 
heading  angles  needed  for  the  computation  of 
the  control  law  are  given  by: 


r  = 


(xx  +  yy)z  -  (*~  ±  ypk 

Vi2+i.-’(i-’+y2  +  z-’)' 
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z  = 


(xy-xy) 
(x2  +  y2 ) 


(25) 
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The  control  chains  defined  in  Expressions  (20) 
and  (21),  together  with  the  expressions  for 
command  generation  and  the  expressions  for 
flight  path  angle  and  heading  angle  are  used  in 
conjunction  with  the  Nonlinear  Synthesis  Tools 
[16]  software  package  to  derive  the  integrated 
flight  control  system.  The  performance  of  this 
integrated  guidance-autopilot  is  then  assessed 
for  an  engagement  scenario  defined  in  the 
following  subsection. 

3.2.1  Engagement  Scenario  1 

In  this  engagement  scenario,  both  the 
warhead  and  target  are  initially  at  an  altitude  of 
45,000  ft,  and  50,000  ft  apart.  Both  vehicles 
have  15-degree  flight  path  angles  and  are  on 
reciprocal  headings,  and  both  have  an  initial 
velocity  of  6000  ft/sec.  In  addition,  the  target 
has  an  initial  offset  of  1000  ft.  in  the  east 
direction. 

The  closed-loop  poles  were  chosen  to  be 
{-99,  -90,  -35,  -30,  -6,  -1}  in  both  channels. 
Trajectories  are  shown  in  Figure  23  and  Figure 
24.  The  terminal  miss  distance  was  0.027  ft. 
The  line-of-sight  rates  are  shown  in  Figure  25. 
It  may  be  observed  that  the  line-of-sight  rate 
corrections  are  a  little  more  distributed  over  the 
engagement  than  in  the  integrated  flight  control 
system  discussed  in  Section  3.2.1.  The  angle  of 
attack  and  angle  of  sideslip  are  shown  in  Figure 
26.  There  is  a  large  initial  acceleration  seen  in 
Figure  27  in  the  lateral  channel  to  eliminate  the 
offset.  The  longitudinal  channel  maintains  some 
acceleration  throughout  the  flight.  A  noticeable 
roll  rate  is  induced  which  can  be  seen  in  the 
angular  velocities  shown  in  Figure  28.  The 
actuator  mass  displacements  and  the  control 
forces  are  shown  in  Figure  29.  As  in  the 
previous  designs,  the  force  requirements  are 
well  within  the  capabilities  of  commercial  off- 
the-shelf  linear  electric  motors. 


Figure  23.  Warhead  and  Target  Trajectories, 
Vertical  Plane 


Figure  24.  Warhead  and  Target  Trajectories, 
Horizontal  Plane 


Figure  25.  Line  Of  Sight  Rates  vs.  Time 
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Time  (sec) 


Figure  26.  Angles  of  Attack  and  Sideslip  vs. 
Time 


Time  (sec) 

Figure  27.  Body-Frame  Inertial  Accelerations 
vs.  Time 


Figure  28.  Body  Angular  Velocities  vs.  Time 
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Figure  29.  Actuator  Mass  Displacements 
vs.  Time 

The  engagement  results  given  in  these 
sections  demonstrate  the  feasibility  of  using  the 
moving  mass  actuation  concept  for  the  kinetic 
warhead.  Note  that  the  miss  distances  obtained 
in  this  preliminary  work  can  be  significantly 
improved  by  introducing  additional  dynamic 
compensation  in  the  guidance-autopilot  loop.  It 
is  important  to  state  here  that  the  present  results 
were  obtained  without  including  various  error 
sources  in  the  guidance-autopilot  loop.  Future 
work  will  examine  the  performance  of  the 
integrated  guidance-autopilot  system  in  the 
presence  of  disturbances  and  sensor  errors,  and 
will  investigate  approaches  for  making  the 
system  robust  with  respect  to  target  maneuvers 
and  other  disturbances.  The  present  research 
shows  that  integrated  guidance-autopilot  system 
design  methodology  is  an  effective  and  flexible 
flight  control  system  design  method. 

A  survey  of  electric  actuators  were 
conducted  during  the  present  research. 
Preliminary  analysis  indicates  that  the  moving- 
mass  control  concept  can  be  realized  using 
readily  available,  commercial,  off-the-shelf 
linear  electric  motor  technology. 

4.  Conclusions  and  Future  Work 

This  paper  described  the  modeling, 
simulation,  and  integrated  guidance  and  control 
of  a  kinetic  warhead  utilizing  moving-mass 
actuators.  Moving  mass  actuation  methodology 
can  be  used  in  any  speed  range  both  in  the 
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atmosphere  as  well  as  outside  it,  as  long  as  there 
is  a  force,  either  aerodynamic  or  propulsive, 
acting  on  the  vehicle.  Since  they  are  contained 
entirely  within  the  airframe  geometric  envelope, 
and  because  no  mass  expulsion  is  involved, 
moving-mass  actuation  technique  offers  several 
advantages  over  conventional  aerodynamic 
control  surfaces  and  reaction  control  systems. 
The  disadvantages  of  the  moving  mass  actuation 
technology  are  that  the  airframe  must  provide 
adequate  internal  space  for  the  moving  masses, 
and  must  have  near-neutral  aerodynamic  static 
stability  characteristics. 

The  present  work  has  developed  a  high 
fidelity,  9  degree-of-freedom  simulation  model 
of  a  kinetic  warhead  with  three  moving-mass 
actuators.  This  simulation  model  was  used  for 
actuator  sizing  and  in  the  development  of  flight 
control  systems.  Commercial,  off-the-shelf 
linear  electric  motors  appear  to  be  capable  of 
providing  the  necessary  force  and  speed  of 
response  characteristics  for  positioning  the 
moving  masses  within  the  airframe. 

A  software  package  for  performing 
numerical  feedback  linearization  technique 
developed  under  a  previous  research  effort  was 
used  to  design  nonlinear  flight  control  systems. 
Pole  placement  technique  was  then  used  to 
design  the  feedback  linearized  flight  control 
system.  Two  different  integrated  design 
methodologies  were  illustrated,  the  first  based 
on  regulating  the  line-of-sight  rates,  and  the 
second  one  controlling  the  trajectory  about  the 
collision  course.  Interception  of  non¬ 
maneuvering  and  weaving  targets  in  both 
atmospheric  and  exo-atmospheric  conditions 
were  illustrated.  In  every  case,  the  moving  mass 
integrated  guidance-autopilot  was  able  to  deliver 
a  miss  distance  smaller  than  half  the  diameter  of 
the  kinetic  warhead.  Additional  improvements 
in  the  guidance-autopilot  system  design  are 
feasible  through  the  introduction  of  dynamic 
compensators  in  the  flight  control  loop. 

Future  research  will  examine  the  impact  of 
uncertainties  in  the  target  maneuvers  and 
environmental  factors  on  the  kinetic  warhead 
performance.  Application  of  the  moving  mass 
control  technology  to  other  flight  vehicles  and 


underwater  vehicles  will  also  be  of  future 
interest. 
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